Field— Controlled Magnetic Order in the Quantum Spin— Ladder System (Hpip)2CuBr4 
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Neutron diffraction is used to investigate the field-induced, antiferromagnetically ordered state 
in the two-leg spin-ladder material (Hpip)2CuBr4. This "classical" phase, a consequence of weak 
interladder coupling, is nevertheless highly unconventional: its properties are influenced strongly by 
the spin Luttinger-liquid state of the ladder subunits. We determine directly the order parameter 
(transverse magnetization), the ordering temperature, the spin structure, and the critical exponents 
around the transition. We introduce a minimal, microscopic model for the interladder coupling and 
calculate the quantum fluctuation corrections to the mean-field interaction. 

PACS numbers: 75.10.Jm; 75.30.Kz; 75.25.+z; 75.40.Mg 



Low-dimensional magnets have been the subject of 
intense theoretical research for many decades. Of par- 
ticular interest are the intriguing ground- and excited- 
state properties of one-dimensional (ID) systems such 
as chains and ladders [TJ|2J|31l31[51[ni[7]. In this con- 
text, residual interactions between the low-dimensional 
units, which are always present in real materials, may be 
viewed as a distraction from the intrinsic physics. How- 
ever, such interactions open up fascinating new avenues 
of investigation concerning the crossover from one- to 
higher-dimensional behavior. In 3D, antiferromagnetic 
(AF) magnons in a gapped, quantum magnet undergo 
Bose-Einstein Condensation (BEG) at a magnetic field 
Be, where the gap is closed by the Zeeman effect [8]. At 
this quantum critical point, the spin components develop 
long-ranged order perpendicular to the magnetic field 
(3D-XY type). By contrast, in ID any long-range order 
is destroyed by quantum phase fluctuations and a critical 
phase with algebraic spin correlations - a spin Luttinger 
liquid (LL) - is predicted ^ [7j . While the spin LL may 
be realized at finite temperatures in coupled S = \/2 
chain systems such as KGuFa [HI HH] , a particularly rich 
phase diagram is expected for weakly coupled ladders: 
here the spin LL is induced from a gapped, quantum dis- 
ordered (QD) phase by an applied magnetic field, and its 
LL parameters can be tuned directly by the field ^Tj. 

Materials realizing quasi-lD spin-ladder geometries. 
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FIG. 1: Low-T phase diagram of (Hpip)2CuBr4. The 
crossover temperature to the spin LL phase is derived from 
MCE measurements and the phase transition to the BEC 
(3D-XY magnetic order) from neutron diffraction. The con- 
tour plot is based on 18 individual field scans of the MCE (two 
shown as gray lines), using {SQ /SB)/T — —{dmz/dT)\B |15) . 
The red line is based on a theoretical fit (see text). 



and with critical fields Be (QD to LL) and Bg (mag- 
netic saturation) accessible in the laboratory, are rare [3]. 
While the latter difficulty is generally overcome in metal- 
organic compounds, these can still suffer from other com- 
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plications, such as the additional terms found in the 
magnetic Hamiltonian of CuHpCl [11] . Thus the system 
(C5Hi2N)2CuBr4 ((Hpip)2CuBr4) [12l [13] is unique in 
its class: numerous thermodynamic measurements are in 
quantitative agreement with predictions for an ideal lad- 
der [in [T3J [T3], and inelastic neutron scattering (INS) 
demonstrates a high degree of one-dimensionality |17j . 
All of these techniques point to a minimal Heisenberg lad- 
der Hamiltonian with respective rung and leg exchange 
constants Jr = 12.9(2) K and Jj = 3.3(3) K. Only nuclear 
magnetic resonance (NMR) measurements have to date 
been performed at temperatures low enough to access the 
energy scale of the inter ladder coupling: in Ref. 16], the 
transition to 3D order is found at and below 110 mK, and 
the phase boundary interpreted in terms of a field-tuned 
spin LL regime between Be and Bg. 

In this Letter we report the results of comprehen- 
sive neutron diffraction studies, performed at dilution- 
refrigerator temperatures, of the 3D ordered phase in 
(Hpip)2CuBr4. We determine the spin structure and 
measure directly the transverse magnetic moment as a 
function of field and temperature. The vanishing of the 
order parameter, combined with the magnetocaloric ef- 
fect (MCE), give the LL crossover, the LL exponents 
close to the 3D regime, and independent measurements 
of the phase boundary. From the spin structure we de- 
duce a microscopic model for the interladder coupling 
which allows a quantitative determination of the inter- 
action between the ID subunits, using a combination of 
mean-field calculations [T5] and Quantum Monte Carlo 
(QMC) simulations of the full interacting-ladder model. 

The MCE was measured on high-quality single- 
crystalline (Hpip)2CuBr4 in a standard dilution refrig- 
erator at the NHMFL in Los Alamos, with sweep 
rates between 0.025 T/min and 0.075 T/min. Neutron 
diffraction experiments were performed on deuterated 
(Dpip)2CuBr4 single crystals with sample mass 200 mg 
on the instruments D23 at the ILL and RITA-2 at SINQ 
(PSI), using standard set-ups. For all measurements, a 
vertical magnetic field was applied along the crystallo- 
graphic 6-axis, i.e. perpendicular to the ladders. 

Before studying the 3D ordered phase, it is necessary to 
understand in full detail the disordered phase from which 
it emerges. The quasi- ID regime above the 3D phase 
boundary is investigated using the MCE, which maps the 
crossover from the QD state into the spin LL through lo- 
cal extrema in the temperature-dependence of the longi- 
tudinal magnetization m^. Figure 1 shows duriz/dT (con- 
tour plot) and these extrema (black circles) down to 100 
mK, extending (from 300 mK) the results of our previous 
measurements |15j and providing a frame of reference for 
other reported results '16]. The LL crossover is analyzed 
by a sliding- window technique ^18j , whose first step is the 
determination of the critical fields Be = 6.96(2) T and 
Ba = 13.85(3) T. With these values fixed, the crossover 
temperature is fitted to Tll (x {B — Bc)'^ (black lines). 
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FIG. 2: Summary of neutron diffraction data, (a) Q-scans 
across an AF Bragg peak after subtraction of a flat back- 
ground measured in the QD phase at B = 6 T and T = 63 
mK. (b) T-dependence of the Bragg intensity, demonstrating 
the onset of 3D long-range order at Tj-i{B); solid lines are fits 
using the 3D-XY exponent, (c) B-dependence of m^, mea- 
sured at g =(1.5 0) for T = 54 mK (blue) and T = 75 mK 
(red). Solid and dashed lines are from DMRG MFA for the 
given interladder interaction J'. Error bars on data points are 
based on counting statistics, while the vertical black line indi- 
cates the systematic error of the calibration to absolute units, 
(d) Magnetic signal at Q =(1 2), which is proportional to 
nil: shown is the neutron intensity after subtraction of the 
nuclear contribution, which also corrects for magnetostriction 
effects. The red line is obtained from a QMC calculation. 



yielding an exponent t/ = 2.1(1) at Be [H]. A ladder 
spin system in a field B ^ Be can be mapped to a free- 
fermion model (LL exponent K = 1), and hence one 
expects = 1 [5]. However, a bosonization interpreta- 
tion of density-matrix renormalization-group (DMRG) 
calculations for the ladder model [T5] shows that K{B) 
decreases rapidly below 1 as the field is moved away from 
Be and B^. Consequently, the true critical regime is very 
narrow, a result seen also in QMC calculations for the 
related Haldane spin chains [50]. Thus while our MCE 
measurements demonstrate spin LL behavior down to 100 
mK, the universal exponents are not reached before the 
LL regime is cut off by 3D order. Close to the critical 
fields, our results indicate that the ladder system remains 
far from this universal regime even for T/Ji w 0.03. 

Turning now to the ordered phase. Fig. 2 summa- 
rizes the results of our neutron diffraction measurements, 
which were taken at temperatures down to 54 mK. Fig- 
ure 2(a) shows Q-scans across the AF wave vector, 
Q =(1.5 0), on cooling the sample from the spin LL 
regime dX Be < B < B^. Resolution-limited magnetic 
Bragg peaks are observed at base temperature, demon- 
strating long-range AF order. The magnetic Bragg peak 
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FIG. 3: Magnetic structure (black arrows) in the 3D ordered 
phase of (Dpip)2CuBr4, determined by neutron diffraction at 
B = 8.6 T and T — 63 mK. Only Cu atoms forming the lad- 
ders are shown (red and blue). The spin component along 
the field (-B||fe) is fixed by QMC calculations. Projections are 
shown (a) on the 6c-plane and (b) on the oc-plane. Interlad- 
der bonds J' for one rung are shown in green. 



remains at the same commensurate position, but its in- 
tensity decreases with increasing field, indicating a sub- 
stantial field-dependence of the transverse ordered mo- 
ment. The temperature-dependence of the Bragg in- 
tensity is presented in Fig. 2(b): the vanishing of the 
magnetic signal is used to determine the phase boundary 
shown as red circles in Fig. 1, and its thermal evolution is 
very well described by the critical exponent, 2/3 = 0.70, of 
the 3D XY model [8 over the full range of data available 
(essentially ^Tn < T < Tn). 

Figure 2(c) shows field scans of the AF Bragg inten- 
sity at Q =(1.5 0) for T = 54 mK and T = 75 mK. 
This is proportional to the square of the transverse mag- 
netization, m^, and was scaled to the ordered moment 
obtained at -B = 8.6 T from a complete refinement of 
the spin structure (below). In contrast to NMR, is 
determined directly by neutron scattering, allowing addi- 
tional quantitative tests of theoretical predictions. From 
Fig. 2(b), the AF order parameter is not saturated at 
T = 54 mK, and no Bragg peak (or 3D transition) is 
observed for _B > 12 T. Above Be, the uniform magne- 
tization TUz increases monotonically, which generates a 
small, ferromagnetic (FM) signal on top of the nuclear 
Bragg peaks. Our characterization of all components of 
the magnetization is completed by field scans of this mag- 
netic intensity, as shown in Fig. 2(d) for Q =(1 2). The 
red line is obtained from a QMC calculation of the ladder 
magnetization, m^, using the exchange interactions cited 
above, and shows good overall agreement. 

In addition to and TOj., neutron diffraction also al- 
lows a quantitative determination of the magnetic struc- 
ture. At base temperature and B = 8.6 T (maximum 



Tn « 110 mK), the intensities of 26 AF Bragg peaks 
were recorded on D23. Among the four allowed mag- 
netic structures, that shown in Fig. 3 provides the best 
fit (x^ = 5.54): the spins are aligned perpendicular to the 
a-axis and antiparallel within the ladder, but parallel on 
ladders of the same type [propagation vector fc = (0.5 
0)]. The ordered moment is 0.41(2) per copper ion. 
Its orientation is parallel to the maximum of the (/-factor 
in the ac-plane [T2]. We note that nearest-neighbor Cu 
atoms between adjacent ladders of opposite type (red and 
blue) have FM ahgned spins [Fig. 3(a)]. 

A qualitative discussion of our results is aided by a spe- 
cific model for the interladder coupling responsible for 3D 
order. Here we restrict our considerations to a minimal 
model with only one interladder interaction parameter, 
of magnitude J' — 0(100 mK). Inspection of the lattice 
structure of (Hpip)2CuBr4, which is similar to that of 
TlCuCls (albeit with the dimer units rather less tilted 
relative to the ladder axis), motivates the choice of the 
bond J3 in the notation of Ref. [H] . This bond connects 
ladders of opposite orientations [red to blue in Fig. 3(a)] 
along the directions ±(1 ±0.5 0.5), meaning that the Cu 
sites are displaced by one unit along the a-axis. AF 
bonds J' = J3 then ensure the FM alignment mentioned 
above, and a completely unfrustrated spin structure. Any 
one ladder rung has four bonds of this type (coordination 
z — 4). INS measurements of the very small triplet dis- 
persions along the b and c axes tend to support this type 
of model; full details of the 3D spin dynamics and ex- 
change paths will be presented elsewhere [17]. 

The measured transition temperature Tn(J5) is shown 
in Fig. 1. The marked asymmetry about B = (Bc+Bg) /2 
{rriz — 1 /2) arises from the changing infiuence of the up- 
per two triplet branches as the field is increased, and 
is in strong contrast to spin systems in two and three 
dimensions, such as BaCuSi206 and TlCuCla [5]. This 
asymmetry is also apparent in QMC studies of coupled 
ladders [22]. The solid red line in Fig. 1 is fitted from 
the theory presented in Ref. [TB]: T^{B) depends on the 
LL parameters of the ladders, which are themselves func- 
tions of B and are obtained from a bosonization interpre- 
tation of DMRG calculations performed for a ladder with 
Jr/Ji = 3.6. Interladder coupling is treated in a mean- 
field approximation (MFA), and the only free parameter 
in Fig. 1 is the magnitude of this exchange constant, J^f- 
The best fit to our experimental data yields J^p = 20(1) 
mK, providing valuable independent confirmation of the 
results of Ref. [16] . 

However, by its nature the MFA underestimates the 
effects of quantum fluctuations, and thus overestimates 
the value of Tn corresponding to a given J' (for coupled 
ID chains by approximately 50% !9,). We have obtained 
the "bare" value of J' by performing QMC simulations 
[211 US] for the interladder exchange geometry of the min- 
imal model, which is a very demanding computational 
task due to the strong spatial anisotropy of the exchange 



4 



- -J'„p= 80 mK 
J' = 50 mK 

MF ^ ^ ^ 
/ 

ll, - - 


— J'^P= aJ'= 
J'„P= aJ'= 
--'''mf=°J- 
s 

^ \ 


59 mK 
37 mK 
20 mK 




III 

1 ^»~^ 


^ ^ 

\ 






■ QMC 


J'=80mK ' 






QMC 


J'= 50 mK 





b °-3 



8 10 12 14 16 
Magnetic field [T] 




100 200 

J' [mKJ 



FIG. 4: The 3D-XY phase boundary by DMRG MFA and 
QMC. (a) Tm(B) for several values of J', as indicated, a = 
0.74(1). (b) Tn as a function of J'. 



constants. It has been shown numerically [23j that quan- 
titative agreement of mean-field and exact results can be 
recovered simply by introducing a renormalization factor 
a, which is a constant in the weak-coupling limit, and 
for the purposes of the current analysis can be viewed 
as a renormalization of the coupling J^p = aJ' . The 
3D ordering temperature computed by QMC is shown in 
Fig. 4(a) at B = 8.9 T (and 11.9 T) for two chosen values 
of J'; the solid lines are obtained from the DMRG MFA. 
As shown in Fig. 4(b), a single value a = 0.74(1) is in- 
deed sufficient to ensure perfect agreement between the 
QMC simulations and the DMRG MFA. Hence the best 
fit to J^p corresponds to a bare (microscopic) coupling 
J' 27(2) mK. 

A further test of the coupling model and the DMRG 
MFA is provided by the transverse-moment measure- 
ments shown in Fig. 2(c). The qualitative shape of the 
zero-temperature DMRG-based results is again asym- 
metric, due to the changing influence of the high-lying 
triplet modes and to the field-dependence of the under- 
lying LL parameters. While this appears to mirror the 
data, a quantitative comparison faces two complications: 
experimentally, the data is not in the low-temperature 
limit, and indeed falls unexpectedly above 10.5 T; theo- 
retically, the calculated moment is rather insensitive to 
the value of J', while the quantum fluctuation suppres- 
sion factor for m^, which need not match that deduced 
from Tn, is unknown. Clearly even the unrenormal- 
ized value, J' = 27 mK, somewhat underestimates the 
transverse moment, and signiflcantly higher values, up 
to J' « 80 mK (Fig. 2(c), P/Tj), appear to be required in 
this framework. More detailed experiments and theoret- 
ical analysis will be necessary to resolve this discrepancy. 

In summary, comprehensive neutron diffraction data 
and measurements of the magnetocaloric effect are pre- 
sented to investigate the 3D ordered phase realized in the 
prototypical spin-ladder material (IIpip)2CuBr4 at low 
temperatures and high magnetic fields. We determine the 
temperature of the transition which separates ID from 
3D physics in coupled spin Luttinger liquids, and charac- 



terize the critical behavior in the ID regime. In the 3D 
phase, we measure both the transverse and longitudinal 
magnetizations, and establish the spin structure. The un- 
conventional field-dependences of the Neel temperature 
and of the transverse magnetization agree well with a 
description based on a minimal coupling model and com- 
bining DMRG calculations with a mean-field treatment 
of the (renormalized) interladder interactions. We deter- 
mine the renormalization factor for this coupled-ladder 
geometry by comparison with detailed QMC simulations. 
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